ABSTRACT The objective was to develop selection indexes for Nellore cattle raised in full-cycle production system in the Brazilian Pantanal. The resulting offspring are retained as replacements or sold at two years of age. Preliminary analyses explored effects of scale on economic values (EV). However, given the available data, these effects were very small. Presented herein are results from a simulated system consisting of 5,000 cows with all animals maintained on pasture as is typical in Pantanal. The EV were determined by approximating the partial derivatives of the profit function, changing one trait at a time, by one unit, while keeping the other traits constant. Traits in the breeding objective were mature cow weight, direct and maternal weaning weight, postweaning average daily gain, subcutaneous fat depth, longissimus muscle area, and stayability. Economic values were calculated on the basis of number of animals (per head), number of animal units, and arroba of carcasss weight. Regardless of the basis, maternal weaning weight and subcutaneous fat depth made negligible contributions to the breeding objective. Proportions of variation in the breeding objectives (per head, per animal unit, per arroba) explained by cow weight, direct weaning weight, postweaning average daily gain, stayability, and longissimus muscle area were: 13, 13, 17; 6, 1, 5; 3, 3, 4; 67, 67, 61; and 11, 17, 13, respectively. These indexes may aid Nellore breeders in their selection decisions, thus facilitating the genetic progress and increased productivity and profitability of Pantanal herds.

Resumo: O objetivo deste trabalho foi avaliar o potencial de diferentes modelos de limiar para a determinação da variabilidade genética em bovinos Nelore, com base nas estimativas de herdabilidade, quanto às características habilidade de permanência no rebanho (STA) e probabilidade de parto aos 36 meses de idade (CP36). Utilizaram-se dados provenientes dos rebanhos da raça Nelore participantes do Programa de Melhoramento Genético Animal da parceria Embrapa-Geneplus. Modelos binomiais e multilimiares foram definidos para as características STA e CP36. As estimativas de herdabilidade foram obtidas com uso de procedimentos bayesianos pelo programa Multiple-trait Gibbs Sampler for Animal Models (MTGSAM), por meio do modelo touro-avô materno. As estimativas de herdabilidade, proporcionadas pelos modelos binário e alternativo, foram, respectivamente, 0,08 e 0,12 para STA e 0,17 e 0,12 para CP36. O modelo de múltiplos limiares é eficiente para a detecção de variabilidade genética quanto à característica habilidade de permanência no rebanho, mas não quanto à probabilidade de parto em vacas com 36 meses de idade.

Abstract: The objective of this work was to evaluate the potential of different threshold models to determine the genetic variability in Nellore cattle, with basis on the heritability estimates for the traits stayability (STA) and first calving probability at 36 months of age (CP36). Data came from the Nellore herds participating in the animal breeding program of the Embrapa-Geneplus partnership. Binomial and multi-threshold models were defined for the STA and CP36 traits. Heritability estimates were obtained following Bayesian procedures in the Multiple-trait Gibbs Sampler for Animal Models (MTGSAM) software, using a sire-maternal grandsire model. The heritability estimates, provided by the binary and alternative models, were, respectively, 0.08 and 0.12 for STA and 0.17 and 0.12 for CP36. The multi-threshold model can efficiently detect the genetic variability for stayability, but not for probability of calving for 36-month-old cows.

Resumo A região Centro-Oeste tem se destacado como produtora de carne bovina e, para tanto, produzir com eficiência tem sido a meta dos produtores. Nesse sentido, objetivou-se estimar parâmetros e tendências genéticas para pesos aos 240 (P240) e 420 (P420) dias de idade em bovinos da raça Nelore, criados na região de Goiás, Brasil. Utilizaram-se 48.580 dados para P240 e 28.685 para P420. Os parâmetros genéticos foram obtidos por análises univariadas em que o modelo continha os efeitos fixos de grupo de contemporâneos (estação de nascimento - águas e seca -, ano de nascimento, sexo e fazenda) e a covariável idade da vaca. Como efeito aleatório, utilizou-se o efeito aditivo direto e materno para P240 e apenas direto para P420. As tendências genéticas foram obtidas por meio da análise de regressão do valor genético do animal sobre ano de nascimento. Os resultados obtidos de herdabilidade direta foram 0,15 ± 0,02 para P240 e 0,24 ± 0,00 para P420. O ganho genético anual foi de 0,274 kg para P240 e de 0,506 kg P420. Os resultados indicam a possibilidade de se obterem maiores ganhos genéticos por meio da seleção, que está sendo mais eficiente sobre o P420.

Abstract The Midwest region has excelled as a beef producer and producing efficiently has been the goal of livestock farmers. This study aimed to estimate the genetic parameters and trends for weights at 240(W240) and 420(W420) days of age in Nellore cattle, reared in Goias state, Brazil. A total of 48,580 data for W240 and 28,685 for W420 were used. Genetic parameters were obtained by univariate analysis. The model included the fixed effects of contemporary group (birth season-dry and wet-, year of birth, sex and farm) and covariate age of cow. Direct and maternal effects were used as random effect for W240 and the direct effect for W420. The genetic trends were obtained by regression analysis of the genetic value of the animal on the year of birth. The results of direct heritability were 0.15 ± 0.02 for W240 and 0.24 ± 0.00 for W420. The annual genetic gain was 0.311 kg for W240 and 0.511 kg for 420 days of age (W420). The results indicate the possibility of obtaining higher genetic gains through selection; however, the selection is being more efficient at W420.

Estudou-se a divergência genética entre 73 reprodutores da raça Gir, por meio de medidas de dissimilaridade, métodos de agrupamento e análises gráficas por componentes principais, com base nas diferenças esperadas em suas progênies relativas a caracteres ponderais e reprodutivos. As distâncias Euclidianas médias padronizadas obtidas foram 4,3743 entre os mais dissimilares e 0,1135 para os mais similares. Quatro grupos de reprodutores com a mesma similaridade foram obtidos por métodos de agrupamento por otimização. Os três primeiros componentes principais explicaram 82,01% da variância total, o primeiro explicou 37,09%, o segundo respondeu por 26,58% e o terceiro por 18,34% da variância. Os escores dos componentes principais possibilitaram a avaliação visual da divergência genética, por meio de gráficos de dispersão. Os resultados encontrados permitiram identificar entre os animais avaliados, os mais divergentes, recomendando que suas progênies possam ser utilizadas como genitores de acasalamentos em programas de melhoramento que visem obter bezerros com melhor desempenho que os pais.

We studied the genetic divergence among 73 Gyr bulls, through dissimilarity measures, grouping methods and graphic analysis by main components based on the differences expected in its progenies related to productive and reproductive characteristics. The mean Euclidian standardized distances obtained were 4.3743 among the most dissimilar ones and 0.1135 for the most similar ones. Four groups of reproducers with the same similarity were obtained through grouping methods by optimization. The three main components explained 82.01% of the total variance; the first one explained 37.09%, the second one responded by 26.58% and the third one by 18.34% of variance. Scores of the main components enabled the visual evaluation of genetic divergence, through dispersion graphics. The results allowed the identification, among the evaluated animals, of the most divergent ones, recommending that their progenies can be used as mating genitors in improvement programs which aim at producing calves with better performance than their parents.

To provide information that can contribute to the genetic evaluation and selection in Tabapuã cattle, genetic parameters were estimated for growth traits using a multi-trait animal model. (Co)variance components were estimated by a Bayesian approach using the Gibbs Sampler. Moderate and similar responses to selection are expected when selecting for weights at the four ages evaluated, since the direct heritability estimates were similar and of average magnitude (0.18 to 0.19). The direct and maternal additive genetic correlations between all pairs of weights were higher than 0.70, indicating a high degree of association between the four traits. This suggests that using any of them as selection criteria will result in correlated response in the others and that multi-trait analysis are recommended for the genetic evaluation of growth traits in beef cattle for exploiting this association.

Objetivou-se estimar parâmetros, correlações e tendências genéticas para os pesos de animais da raça Nelore criados na região da Mata e no Agreste Nordestino aos 205 (P205), 365 (P365) e 550 (P550) dias de idade. Para realização das análises genéticas, utilizou-se da metodologia de modelos mistos. O modelo estatístico continha os efeitos aleatórios aditivos diretos, aditivos maternos, ambiente permanente e o erro; efeitos fixos de grupo contemporâneo - fazenda, sexo, estação (água e seca), mês e ano de nascimento do bezerro e a covariável idade da vaca ao parto. As tendências genéticas foram estimadas pela regressão dos valores genéticos sobre o ano de nascimento dos animais. As estimativas da herdabilidade direta foram 0,13 ± 0,03 para P205; 0,03 ± 0,02 para P365 e 0,12 ± 0,04 para P550. As correlações variaram de 0,55 a 0,85. Os parâmetros genéticos permitem ganhos por meio da seleção e, de acordo com as correlações, a seleção para peso a determinada idade possibilitará ganhos as idades posteriores, podendo assim realizar-se seleções em idades mais jovens, apesar dos ganhos serem lentos devido aos baixos coeficientes de herdabilidade.

The study aimed to was estimate the genetics parameters, correlations and trends for the weight of Nellore breed bovines at 205, 365 and 550 days of age. Genetics analysis was carried out by the methodology of mixing models, the statistical model containing the aleatory direct and maternal genetic effects, the uncorrelated maternal permanent environmental effect and the error, as fixed effects of contemporary group (farm, sex, season, water and drought, the year of birth of the calf), and the covariable age of dam. The genetic trends were estimated by linear regression of the breeding values on years of birth of the animals. Estimates of direct heritabilities were 0.13 ± 0.03 to W205; 0.03 ± 0.02 to W365; and 0.12 ± 0.04 to W550. Genetic correlations were of 0.55 a 0.85. The genetic parameters allow gains through selection and, according to the correlations, selection for weight gain to a certain age will allow gains weight at later ages, despite the gains are slow due to low heritability coefficients.

The objective of the present study was to analyze models with genetic and/or residual heteroskedasticity for genetic evaluation of the weaning weight of Charolais-Zebu crossbred calves. Weaning weight data from 56,965 crossbred calves were analyzed using animal models with different combinations of genetic and residual heteroskedasticity and/or homoskedasticity. The inference on a posteriori distributions of genetic parameters were by the Monte Carlo method via Markov chains. The model with genetic and residual heteroskedasticity was the best fit on the data. Groups of animals with different genetic compositions, expressed as percentages of Charolais-Zebu breed alleles and individual and maternal heterozygosis, had different genetic variances. These genetic variances could be modeled by linear functions of the Charolais and Zebu genetic variances and the variance attributed to segregation. The breed compositions, the individual and maternal heterozygosis, the sex and age of dam at calving were significant sources of residual heteroskedasticity. The a posteriori means for heritabilities and sire and dam classifications were altered due to genetic and residual heteroskedasticity.

Este trabalho foi realizado com o objetivo de avaliar alternativas para modelagem de efeitos genéticos e ambientais que influenciam o peso à desmama de bezerros mestiços Charolês-Zebu. Foram utilizados os dados de peso à desmama de 56.965 bezerros, com percentuais de alelos de origem da raça Charolesa variando entre 23% e 66%. Os modelos considerados diferiram quanto à inclusão de efeitos genéticos aditivos e não-aditivos médios, individuais e maternos, bem como quanto à maneira de modelar estes efeitos (regressão múltipla ou variável classificatória). Também foi avaliada a inclusão da interação entre o grupo genético do bezerro e seu grupo de contemporâneos, como um efeito aleatório não correlacionado. De acordo com os critérios de ajuste utilizados, a interação entre o grupo genético e o grupo de contemporâneos representa parte significativa da variação observada para o peso à desmama de bezerros Charolês-Zebu e os efeitos dos percentuais de alelos de origem da raça Charolesa e de heterozigoses, individuais e maternos, podem ser modelados por meio da utilização de regressão múltipla.

This study was carried out to evaluate some alternatives for modeling of genetic and environmental effects influencing the weaning weight of Charolais-Zebu crossbred calves. The weaning weight data of 56,965 calves with alleles' percentage from the Charolais breed varying from 23% to 66% were used. The models under comparison differed for the inclusion of the individual and maternal additive and non-additive genetic effects, as well as the way to model those effects (multiple regression or discrete variable). The inclusion of the interaction between the calf genetic group and its contemporary group was also evaluated as a random uncorrelated effect. The adjustment criteria pointed out the interaction between the calf genetic group and its contemporary group to represent a significant part of the variation observed for the weaning weight of the Charolais-Zebu crossbred calves, and the effects of alleles' percentage from the Charolais breed and the heterozygosity, individual and maternal, can be modeled by using the multiple regression.

O objetivo deste trabalho foi estimar parâmetros genéticos para idades ao primeiro (IPP) e segundo partos (ISP), idades ao primeiro (IPP PEN) e segundo partos penalizadas (ISP PEN), e a ocorrência de parto até os 38 meses de idade (PP38), em bovinos da raça Canchim, bem como estimar as relações genéticas entre essas características e o peso (PE420) de machos e femêas e a circunferência escrotal (CE420) de novilhos, medidos aos 420 dias de idade. Foram estimados parâmetros genéticos em análises bicaracterísticas, por meio de inferência bayesiana. As análises foram realizadas entre CE420 e PE420 e as características reprodutivas medidas em fêmeas. As características reprodutivas das fêmeas devem apresentar baixa resposta à seleção, e mudanças no manejo e nas condições ambientais melhoram tais características. A seleção quanto ao maior peso aos 420 dias de idade não deve provocar mudanças nas características reprodutivas em fêmeas, enquanto a seleção para o aumento da circunferência escrotal, na mesma idade, contribui para melhorar tais características.

The objective of this work was to estimate genetic parameters for age at first (AFC) and at second (ASC) calvings, penalized age at first (AFC PEN) and at second (ASC PEN) calvings, and occurrence of calving until 38 months of age (PP38) in Canchim breed, as well as to estimate the genetic correlations between these traits with weight (PE420) of males and females, and scrotal circumference (CE420) measured at 420 days of age. Genetic parameters were estimated by two-trait analyses using the Bayesian inference. The analyses were performed between CE420 and PE420, and the reproductive traits measured in the females. The female reproductive traits should show low response to selection, and changes in the management and environmental conditions improve these traits. Selection for the higher weight at 420 days of age does not change female reproductive traits, while selection to increase scrotal circumference, at the same age, contributes to improve such traits.

Os efeitos de ano e época de nascimento, propriedade, sexo do bezerro, regime alimentar, idade do bezerro e idade da vaca ao parto, os efeitos aditivos direto e materno da raça Charolesa e o efeito de heterose materna entre as raças Charolesa e zebu foram avaliados pelo método dos quadrados míninos sobre características de importância econômica, como peso, perímetro escrotal e escores visuais de conformação frigorífica, de umbigo e de pelame à desmama em bovinos da raça Canchim. Utilizaram-se dados de 12.334 animais Canchim nascidos de 1999 a 2005. O ano e a época de nascimento, a propriedade, o regime alimentar e a idade da vaca ao parto (linear e quadrático) tiveram efeitos significativos sobre todas as características estudadas. Houve efeitos do sexo sobre o peso do bezerro à desmama e o escore de umbigo; da idade do bezerro sobre peso, perímetro escrotal e escores de conformação frigorífica e de pelame; e das proporções de Charolês no bezerro e na mãe sobre o peso, a conformação frigorífica e a qualidade de pelagem; e da heterozigose materna sobre peso, perímetro escrotal e escores de conformação frigorífica e de umbigo. Os efeitos aditivos diretos da raça Charolesa como desvio do zebu foram negativos sobre peso, conformação frigorífica e qualidade da pelagem à demama, enquanto os efeitos aditivos maternos foram positivos. Os efeitos heteróticos maternos foram positivos sobre peso, perímetro escrotal e escores de conformação frigorífica e de umbigo. Em geral, essas fontes de variação devem ser consideradas nas estimativas de parâmetros genéticos e nas avaliações genéticas dos animais Canchim para as características estudadas.

The effects of calving year and season, property, sex of the calf, feeding regime, age of calf and age of dam at calving, maternal and direct additive effects of Charolais breed, and maternal heterosis effects among Charolais and zebu breeds were evaluated by the method of least squares on economical importance traits such as weight, scrotal circumference and visual scores for slaughter conformation, navel score and hair coat at weaning of the Canchim breed. Data from 12,334 Canchim calves born from 1999 to 2005 were used. Calving year and season, property, feeding regime and age of dam (linear and quadratic) had significant effects on all the studied characteristics. There were effects of sex on weaning weight and navel score; of age of the calf on weight, scrotal circumference and slaughter conformation and hair coat scores; and of proportion of Charolais on the calf and on the dam on weaning weight, and on scores of slaughter conformation and hair coat; and of maternal heterozygose on weight, scrotal circumference, slaughter conformation and navel scores. Additive direct effects of the Charolais as deviation of zebu were negative on weight, slaughter conformation and hair coat score at weaning, while the maternal additive effects were positive. Maternal heterotic effects were positive on weight, scrotal circumference, slaughter conformation and navel scores. In general, these sources of variation have to be taken into account on estimates of genetic parameters and on genetic evaluations of Canchim animals for the studied characteristics.

O objetivo deste trabalho foi estabelecer o modelo mais adequado para avaliação genética de bovinos Canchim e estimar os parâmetros genéticos de características produtivas à desmama. Foram utilizados dados de: 12.103 animais, quanto ao peso (PD); 5.278, quanto ao perímetro escrotal (PE); 8.343, quanto ao escore visual da conformação frigorífica (CF); 9.111, quanto ao escore de umbigo (UM); e 7.986, quanto ao escore de pelame (PEL). Os modelos estatísticos incluíram os efeitos fixos e os efeitos aleatórios genéticos aditivos direto, materno e de ambiente permanente materno, em diferentes combinações. As análises foram feitas pelo método da máxima verossimilhança restrita livre de derivadas. O modelo completo foi o mais adequado para PD, PE, CF e UM, enquanto o modelo com uso apenas dos efeitos genéticos aditivos direto e materno foi o mais adequado para PEL. As estimativas de herdabilidade direta foram 0,17, 0,13, 0,20, 0,18, e 0,52 para PD, PE, CF, UM e PEL, respectivamente, o que indica a possibilidade de se obter progresso genético por meio da seleção para essas características, principalmente para PEL. As correlações genéticas aditivas diretas entre as características variaram de -0,16 a 0,61. As correlações entre PD e PE e entre PD e CF indicam que a seleção para PD deve proporcionar ganho genético em PE e CF

The objective of this study was to establish the best-fit model for the genetic evaluation of Canchim cattle and to estimate the genetic parameters for productive characters at weaning. The data used were related to: 12,103 animals, for weight (PD); 5,278, for scrotal circumference (PE), 8,343, for slaughter conformation visual score (CF); 9,111, for sheath and navel scores (UM), and 7,986, for hair coat score (PEL). Statistical models included fixed effects and random additive direct, maternal and maternal permanent environmental effects, in different combinations. The analyses were done by the derivative-free restricted maximum likelihood method. The complete model was the most suited for PD, PE, CF and UM, while the model with only additive direct and additive maternal random effects was the best for PEL. Direct heritability estimates were 0.17, 0.13, 0.20, 0.18, and 0.52 for PD, PE, CF, UM and PEL, respectively, which indicates that it is possible to obtain genetic progress through selection for these traits, mainly for PEL. The additive direct genetic correlations among traits varied from -0.16 to 0.61. The correlations between PD and PE and between PD and CF indicate that selection for PD should result in genetic gain in PE and CF

Este trabalho foi realizado com o objetivo de avaliar alternativas para modelar a interação entre idade ao parto e grupo genético da vaca (GGV) sobre o peso à desmama de bezerros mestiços Charolês-Zebu. Dados de peso à desmama de 56.965 bezerros Charolês-Zebu foram analisados utilizando o método de quadrados mínimos. Foram estimados coeficientes de regressão para idade da vaca ao parto (IVP) dentro de cada grupo genético da vaca (modelo CLA); para IVP × fração esperada de alelos de origem da raça Charolesa na vaca (IVP × FCh) e IVP × fração esperada dos loci da vaca com um alelo proveniente da raça Charolesa e outro proveniente de raças zebuínas (IVP × FH) (modelo FChFH); para IVP × FCh (modelo FCh); para IVP × FH (modelo FH); ou apenas para IVP (modelo SINT). A forma geral do efeito da IVP sobre o peso à desmama foi modelada por polinômios segmentados com efeitos linear e quadrático gerais e quadráticos específicos, a partir de cada nó (6,33 e 10,66 anos). Os coeficientes de regressão foram estimados dentro de cada sexo do bezerro em todas as situações. De acordo com o teste F para redução da soma de quadrados do resíduo, a inclusão da interação IVP × FH não melhorou significativamente os ajustes e o modelo CLA foi o que melhor se ajustou aos dados. Contudo, para classes de GGV × sexo do bezerro com reduzido número de observações, o modelo CLA proporcionou estimativas inconsistentes para o efeito da interação IVP × GGV e o modelo FCh mostrou-se adequado. Para as classes de GGV × sexo do bezerro com grande número de observações, as diferenças entre as estimativas do efeito da interação IVP × GGV obtidas pelos modelos CLA e FCh foram de pequena magnitude.

The objective of this study was to evaluate alternatives for modeling the interaction between age of dam at calving (AOD) and the dam genetic group (DGG) on the weaning weight (W225) of Charolais-Zebu (Ch-Z) crossbred calves. Data from 56,965 crossbred calves were analyzed by the least square method. Regression coefficients for age of dam at calving were estimated nested into each class of the dam genetic group (CLA model); for age of dam at calving × dam Charolais percentage (age of dam at calving × FCh) and age of dam at calving × dam heterozygosity (age of dam at calving × FH) (FChFH model); for age of dam at calving × dam Charolais percentage (FCh model); for age of dam at calving × FH (FH model); or only for age of dam at calving (NINT model). Segmented polynomials were used to model the general shape of the age of dam at calving effect and its interaction with dam genetic group. The knots were at 6.33 and 10.66 years of age of dam at calving and general linear and quadratic coefficient regression and specific quadratic coefficient regression after each knot were fitted. The regression coefficients were estimated nested within sex of the calf in all situations. According to the F test for sum of squared residuals differences, the inclusion of the age of dam at calving × FH interaction did not improve the fit of the model and the CLA model provided the best fit. However, the estimates of the age of dam at calving and dam genetic group interaction from the CLA model for dam genetic group × sex of the calf classes with few records were not appropriate, but the estimates of the age of dam at calving and dam genetic group interaction from the FCh model for those classes were appropriate. The differences were small in the estimates of the age of dam of calving and dam genetic group interaction from the CLA or FCh models for dam genetic group × sex of the calf classes with many records.

Este trabalho foi realizado com o objetivo de avaliar alternativas para a modelagem do efeito da idade da vaca ao parto sobre o peso à desmama de bezerros mestiços Charolês-Zebu. Os dados de 56.965 bezerros foram analisados por meio de modelos estatísticos, considerando os efeitos fixos de grupo de contemporâneos, grupos genéticos do touro e da vaca e da idade da vaca ao parto. O efeito da idade da vaca ao parto foi modelado utilizando-se classes anuais de idade, polinômios ordinários (PO) de ordem quadrática até quíntica, ou polinômios segmentados com 2, 3, 4, 6 ou 12 intervalos do mesmo tamanho. No caso dos polinômios segmentados, consideraram-se os efeitos linear e quadrático gerais e apenas o quadrático a partir de cada nó. Os efeitos da idade da vaca ao parto foram considerados aninhados em sexo do bezerro em todos os casos. De acordo com os critérios de ajuste utilizados, ou seja, o teste F para diferenças de soma de quadrados, coeficiente de determinação do modelo, soma de quadrados do resíduo e média dos resíduos ao quadrado, o polinômio segmentado com três intervalos (dois nós) se ajustou tão bem aos dados quanto os polinômios mais complexos.

The objective of this study was to evaluate alternatives for modeling the age of dam at calving (AOD) effect on the weaning weight of Charolais-Zebu crossbred calves. Data from 56,965 calves were analyzed, using statistical models considering the fixed effects of the contemporary groups, sire and dam genetic groups, and AOD. The AOD effect was fitted to models using annual age classes, and ordinary quadratic to quintic-ordered polynomials (OP) or segmented polynomials (SP) with two, three, four, six and twelve evenly spaced intervals. In the case of segmented polynomials, general linear and quadratic effects and only one quadratic additional term from each knot were considered. The AOD effects were nested within sex of calf in all cases. According to the fitting criteria, the F-test for the reduction of residual sum of squares, coefficient of determination, residual sum of squares and mean of squared residuals, the three interval segmented polynomial (two knots) fitted to the data as well as the more complex polynomials.

Foram analisados 117.082 registros de lactações encerradas para a produção de leite corrigida para os 305 dias de lactação (PL305) de 49.676 vacas da raça Holandesa, provenientes de 308 rebanhos distribuídos em sete bacias leiteiras no estado do Paraná, com o objetivo de verificar a existência de interação genótipo × ambiente para a PL305 desses animais utilizando-se a inferência Bayesiana. Todos os animais foram controlados oficialmente entre janeiro de 1992 a dezembro de 2003 pelo Serviço de Controle Leiteiro Mensal da Associação Paranaense de Criadores de Bovinos da Raça Holandesa. Os componentes de co-variância e os parâmetros genéticos foram estimados por meio de análises uni e multicaracteres, de modo que, na análise multicaracter, a PL305 em cada uma das bacias foi tratada como uma característica diferente. A produção de leite corrigida para os 305 de lactação, em kg, nas bacias leiteiras de Castro, Carambeí, Witmarsum, Arapoti, Sul, Norte e Oeste foram de 8.414 ± 1.825, 8.481 ± 2.010, 7.636 ± 1.594, 7.850 ± 1.795, 8.617 ± 2.050, 7.401 ± 1.809 e 7.336 ± 2.456, respectivamente. A estimativa de herdabilidade mais alta (0,39) foi obtida para a bacia leiteira do Oeste e a mais baixa (0,23) para a de Carambeí. As correlações genéticas obtidas entre as bacias leiteiras foram baixas (0,09 a 0,57). As correlações de Pearson e de Spearman mais baixas foram obtidas para a bacia leiteira do Oeste do Paraná e variaram de 0,37 a 0,41 e de 0,37 a 0,49, respectivamente. Esses resultados comprovam a existência de interação genótipo × ambiente nas bacias leiteiras do estado do Paraná.

It was analyzed 117,082 complete lactation records adjusted for 305 days in milk yield (MY305) of 49,676 Holsteins cows from 308 herds distributed among seven dairy production units in the state of Paraná, Brazil. The objective was to verify the genotype × environment interaction to MY305 of Holsteins cows between dairy production units in the state of Paraná using the Bayesian inference. All animals were officially controlled by the Serviço de Controle Leiteiro Mensal of Associação Paranaense de Criadores de Bovinos da Raça Holandesa, between January/1992 and December/2003. The (co)variance components and genetic parameters were estimated using one and multiple-trait analysis, where MY305 of each dairy production unit was considered as a different variable. Means and standard deviations for MY305 (kg) for dairy production units of Castro, Carambeí, Witmarsum, Arapoti, Sul, Norte and Oeste, were 8,414 ± 1,825; 8,481 ± 2,010; 7,636 ± 1,594; 7,850 ± 1,795; 8,617 ± 2,050; 7,401 ± 1,809 and 7,336 ± 2,456, respectively. The highest heritability estimated was obtained for the Oeste unit (0.39) while Carambeí unit presented the lowest value (0.23). Genetic correlations were low (0.09 to 0.57) between dairy production units. Pearson and Spearman correlations were lower for the Oeste unit, when compared with other dairy production units, which ranged from 0.37 to 0.41, and from 0.37 to 0.49, respectively. These results show the presence of genotype × environment interaction between dairy production units in the state of Paraná.

Estimaram-se as herdabilidades e as correlações genéticas para as produções de leite (PL), de gordura (PG) e de proteína (PP) e para os percentuais de gordura (%G) e de proteína (%P) no leite de vacas da raça Holandesa no estado do Paraná. Foram analisados 117.082 registros de lactações encerradas para estimativa das características PL, PG e PP e 116.747 registros de lactações para %G e de %P. Todos os animais foram controlados oficialmente pelo Serviço de Controle Leiteiro Mensal da Associação Paranaense de Criadores de Bovinos da Raça Holandesa entre janeiro de 1992 e dezembro de 2003. A estimação dos componentes de co-variância e dos parâmetros genéticos foi realizada por meio de inferência Bayesiana em análises uni e tricaracteres utilizando-se um modelo animal que incluiu os efeitos fixos de grupo de contemporâneos e número de ordenhas, além dos efeitos linear e quadrático da idade da vaca ao parto e os efeitos aleatórios genéticos aditivos, permanente de ambiente e residual. As médias e os respectivos desvios-padrãopara as PL, PG e PP, em kg, foram 8.181,23 ± 1.903,51; 270,88 ± 67,58 e 249,01 ± 55,65, respectivamente, e para os %G e %P foram 3,33 ± 0,45 e 3,06 ± 0,19, respectivamente. As estimativas de herdabilidade para PL, PG, PP, %G e %P obtidas nas análises unicaracteres foram 0,26; 0,28; 0,25; 0,60 e 0,58, respectivamente, e indicam que as características possuem variação genética aditiva que garante boa resposta à seleção. As estimativas das correlações genéticas entre PL e PG, PL e PP, PG e PP, PL e %G, PL e %P, %G e %P foram 0,56; 0,89; 0,67; -0,42; -0,39 e 0,56, respectivamente, portanto, a seleção realizada para aumentar a PL promoverá aumento na PG e PP e redução nos %G e %P.

The objectives of this study were to estimate heritabilities and genetic correlations for milk (MY), fat (FY) and protein (PY) yields and for percentages of fat (%F) and protein (%P) in milk of Holstein cows in Paraná State. Data consisted of 117,082 lactation records for MY, FY and PY and 116,747 records for %F and %P from herds supervised by the Milk Recording Service of the Holstein Breeders Association of Paraná from January 1992 to December 2003. Covariance components and genetic parameters were estimated by Bayesian inference using an animal model including the fixed effects of contemporary group, number of milking, linear and quadratic effects for calving age and the additive genetic, permanent environmental and residual random effects. Averages and respective standard deviations for MY, FY and PY (kg) were respectively 8,181.23 ± 1,903.51, 270.88 ± 67.58 and 249.01 ± 55.65 and for %F and %P they were 3.33 ± 0.45 and 3.06 0.19, respectively. Heritability estimates for MY, FY, PY, %G and %P obtained from single-trait analyses were 0.26, 0.28, 0.25, 0.60 and 0.58, respectively. They indicate that these traits may be improved by selection. Genetic correlation estimates between MY and FY, MY and PY, FY and PY, MY and %F, MY and %P and for %F and %P were 0.56, 0.89, 0.67, -0.42, -0.39 and 0.56, respectively. Selection for increasing MY will result in positive correlated genetic gains for kg of fat and protein but decrease their percentages in milk.